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Abstract. We describe the formation of a seed massive black hole (MBH) inside a supermassive star (SMS) in a 
distant galactic nucleus. The short-lived SMS is naturally formed due to collision destructions of normal stars in 
the evolving galactic nucleus. The neutron stars (NSs) and stellar-mass black holes form a compact self-gravitating 
subsystem deep inside a SMS. This subsystem is short-lived in comparison with a host SMS and collapses finally 
into the MBH. Just before gravitational collapse of compact subsystem the frequent NS collisions are accompanied 
by the generation of numerous ultra-relativistic fireballs. A combined ram pressure of multiple coexisting fireballs 
produces a quasi-stationary rarefied cavity in the central part of SMS. The protons are accelerated in the fireballs 
and by relativistic shocks in the cavity. All secondary particles, produced in collisions, except the high-energy 
neutrinos are absorbed in the SMS interiors. An estimated high-energy neutrino signal from this hidden source 
can be detected by the neutrino telescope with an effective area 5* 1 km'^ providing the evidence for MBH 
formation in a distant galactic nucleus. A corresponding lifetime of this high-energy hidden neutrino source is 
-0.1-1 yr. 
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1. Introduction 

The origin of massive black holes (MBHs) in the galactic 
nuclei is a long standing problem. In this paper we de- 
scribe a 'natural' scenario of a seed MBH formation in the 
process of dynamical evolution of a central stellar cluster 
in the galactic nuclei. We demonstrate that high-energy 
neutrinos may be the observational signature of MBH for- 
mation. 

High-energy (HE) neutrino radiation from astrophys- 
ical sources is accompanied by other types of radia- 
tion, most notably by the HE gamma-radiation. This 
HE gamma-radiation put an upper limit on the neu- 
trino flux emitted by the neutrino-transparent astrophys- 
ical sources. For example, if neutrinos are generated by 
the interaction of HE protons with low energy pho- 
tons in the extragalactic space or in gamma-transparent 
sources, the upper limit on diffuse neutrino flux I^{E) 
can be derived from electro-ma gnetic cascade radiation 
ijBerezinskv and Smirnovl. Il975|) . However, there can be 
the "hidden" neutrino sources where accompanying X- 
ray and gamma-radiations are absorbed. There are sev- 
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eral known models of the hidden HE neutrino sources : 
The Young SN shells ijBerezinskv and Prilutskvl Il987j) . 
opaque for all radiation except of neutrino during a 
time interval ti, ^ 10'^ — 1 0*^ s. The Thome- Zytkow star 
l|Thorne and ZvtkowL Il977|) . which is a neutron star sub- 
merged i nto a red giant. The cocoo ned massive black hole 
in AGN f|Berezinskv and Ginzbur3 . [198 1,1 where the elec- 
tromagnetic radiation is absorbed in a cocoon around the 
black hol e. The AGN with a st anding shock in the vicinity 
of MBH (IStecker et all Il99ll) can produce large flux of 
HE neutrinos accompanied by the relatively small flux of 
X-ray radiation. 

Recently we proposed the model for a powerful but 
short-lived hidden neutrino source which originates in a 
distant galactic nucleus prior to its collapse into the MBH 
ijBerezinskv and Dokuchagvl l200l|) . In this model we as- 
sume that MBH is formed by the natural dynamical evo- 
lution of the central stellar cluster in the normal galactic 
nucleus. Dynamical evolution of dense central stellar clus- 
ters in the galactic nuclei is accompanied by the growth of 
a central gravitational potential (f> ^ v'^, where w is a veloc- 
ity dispersion of constituent stars. This evolution process 
is accompanied by the contraction of a central part of the 
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galactic nucl eus and term inated by formation of MBH (see 
for a review iReed l)l984l) and references therein). During 
this evolution a dense galactic nucleu s inevitably proceeds 
through the stellar collision phase JSpitzer and Saslawl 
Il96fit ISnitzel Il97ll: iReel Il984l: iDokuchaevi Il99lh . when 
normal stars in the cluster are destroyed in mutual colli- 
sions. A direct collision of two stars ends with their com- 
plete destruction if their relative kinetic energy exceeds 
the gravitational bound energy of stars. This condition 
is realized in stellar clusters with a high enough veloc- 
ity dispersion. Namely, v > Wcsc = (2G'TO*/r,)^/^ is ful- 
filled, where Vesc is an escape velocity from the surface of 
star with mass m* and radius r*. For a solar type star 
Uosc — 620 km s~^. In a stellar cluster with v > v^sc 
the normal stars are eventually destroyed in mutual col- 
lisions or in collisions with the extremely compact stellar 
remnants: neutron stars (NSs) or stellar mass black holes. 
Only compact stellar remnants survive through the nor- 
mal star-destruction stage of evolution at w ~ Vesc and 
form the self-gravitating subsystem surrounded by a mas- 
sive gas envelope. We shall refer to all compact stellar rem- 
nants as the NSs for brevity. The fireballs after numerous 
NS collisions in the contracting NS cluster result in the 
formation of a rarefied cavity in the massive gas envelope. 
Particles accelerated in this rarefied cavity interact with 
a dense envelope and produce secondary HE neutrinos. 
The accompanying gamma-radiation is fully absorbed in 
the case of thick enough envelope. The resulting hidden 
HE neutrino source is short-lived and very powerful: neu- 
trino luminosity may exceed the Eddington limit for the 
electromagnetic radiation. 



IBerezinskv and DokuchaevI l)200l|) considered the case 
when a self-gravitating cluster of NSs in the evolving 
galactic nucleus is formed simultaneously with a formation 
of massive gas envelope. It is possible because the time of 
collision of self-destruction for normal stars appears to be 
of the same order as the time of NS cluster formation at a 
star-destruction stage of evolution of the galactic nucleus 
(see Section O for more details). In this case the newly 
formed NS cluster evolves much faster than a surround- 
ing massive envelope. The corresponding parameters of 
massive envelope are fixed at the moment of NS cluster 
formation when v ~ Wosc- Here we will consider another 
and probably more natural case, when a self-destruction 
of normal stars and formation of a massive envelope in 
the evolving galactic nucleus proceeds faster than a for- 
mation of the self-gravitating subsystem of NSs in it. This 
case corresponds to the initial formation of a supermas- 
sive star (SMS) with the individual NSs submerged in it. 
Correspondingly, a formation of self-gravitating subsys- 
tem of NSs is delayed with respect to the formation of 
SMS, and the evolutionary contraction of SMS must be 
taken into account. We show that the hidden source of 
HE neutrino arises in this cas e also and even more power- 
ful th an in the first scenario ijBerezinskv and DokuchaevI 
l200l|l of a simultaneous formation of a massive envelope 
and NS cluster in the evolving galactic nucleus. 



2. Supermassive star formation 

The SMS may be formed from the gas produced in the pro- 
cess of destructive collisions of stars in the evolved stellar 
cluster with a velocity dispersion v > Wcsc- A characteris- 
tic time-scale for stellar cluster dynamical evolution is the 
(two-body) relaxation time 
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4.6 X 10*7V, 



yr, (1) 



where N = lO'^A^s is the number of stars in the cluster, 
A = ln(0.4iV) is the Coulomb logarithm, w is a star veloc- 
ity dispersion, n is a star number density, m ~ Mq is a 
mass of constituent stars. The corresponding virial radius 
of the stellar cluster is i? = GNm/2v^ . The first equality 
in is valid for the local values of parameters, mean- 
while, the second one is valid only for the mean (virial) 
parameters of a self-gravitating cluster. At > Vesci where 
Vesc is an escape velocity from the surface of constituent 
star, the time-scale for self-destruction of no rmal stars in 
mutu al collisions is tco\i = {vesc/v)'^Ati- (e. g. iDok^ichaevl 
Numerical modelling of ca tastrophic stellar coll i- 
sions has been p erformed by e. g. iBenz and Hilil l)l992|) : 
iLai et al.l l)l993|) . We choose characteris- 
tic threshold value for a complete destruction of two stars 
and final production of unbound gas cloud. As a result, the 
stellar cluster in the evolved galactic nucleus with v > Vesc 
transforms finally into the SMS due to catastrophic stellar 



collisions. At v 



the time-scale for the formation of 



SMS due to destructive collisions of stars is of the same 
order as the relaxation time, tcoii(i'esc) ~ ir('i'esc)- A total 
mass of the gas produced by destruction of normal stars 
composes the major part of a progenitor central stellar 
cluster in the galactic nucleus. Thus, the natural range of 
masses for the formed SMS is Msms = 10^ - lO'^M©. 

A newly formed SMS with mass Msms and radius 
RsMS gradually contracts due to radiation with the 
Eddington luminosity Le = 47rGAfsMS'7tpc/(7T, where uip 
is the proton mass and ctt is the Thompson cross-section. 
A nonrotating SMS is a short-lived object that collapses 
due to post-Newtonian instability. Rotation provides the 
stabilization of SMS if the rotation energy is an apprecia- 
ble part of its total energy Esms — {GMg^g/2RsMs)- In 
general an evolution time of SMS is the Kelvin-Helmholtz 
time-scale isMS — Esms/Le ^ ^sms- Stabilization of 
SMS by rotation (and additionally by internal magnetic 
field) ensures in principle its gradual contraction up to 
the gravitation al radius ('Zel'dovich and Novikovi Il97lt 
IShapiro andTeukolskv. 19M; New and Shapira l200lj) . 
Thus, the maximum evolution time of SMS is of the order 
of the Eddington time = 0.1MsMsc^/-f^E — 4.5 x lO"^ yr. 
We approximate the distribution of gas in a SMS by poly- 
tropic model with an adiabatic index 7 = 4/3. For this 
value of adiabatic index the central gas density in SMS is 
Pc — fcc"-SMS™p and central sound velocity Cs,c — kgVsus, 
where nsMS is a SMS mean number density, vsms — 
(GMsMs/2i?sMs)^^^ is a SMS virial velocity and numeri- 
cal constants kc — 54.2 and kg — 1.51 respectively. 
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3. Neutron star cluster formation 

The compact stellar remnants in the form of NSs and stel- 
lar mass black holes submerged into the newly formed 
SMS. This is a specific feature of SMS formed in the evolv- 
ing stellar cluster. Let SMS contains a small fraction of 
identical NSs with a total mass Mns = /ns-^smsj /ns ^ 
1 and a typical mass of individual NS is tons = 1-4Mq. 
In numerical estimations we will use /ns = 10^^/-2 with 
/_2 ^ 1- An individual NS with a mass m and local ve- 
locity V is spiraling down the S MS center under influ- 
ence of the dyn amical friction force l)Chandrasekhail|l943l: 
lOstrikeil K^m Fdf = A7rI{Gm)^ psus/V^ , where psus is 
a SMS gas density. The dimensionless factor here is / ~ 1 
in the case of V ~ Cg, where Cg is a sound velocity in 
SMS. Deep inside SMS the individual NSs are spiralling 
down and combine into a fast evolving self-gravitating 
subsystem (NS cluster). A corresponding time-scale for 
the NS dynamical friction drag toward the SMS center 
is tdf = V/V = mV/Fdi. At the moment of SMS for- 
mation, when V ^ Cs ^ fosc, a time-scale of NS friction 
drag is of the same order as the time-scale of SMS for- 
mation, tdt{vesc) ~ ^coll(^'esc) ~ ir(wosc)- Thus, in gen- 
eral it is possible either (i) a simultaneous formation of 
a self-gravitating cluster of NSs and the host SMS or 
(ii) a delayed formation of a self-gravitating cluster of 
NSs inside the host SMS. The ca se (i) was considered by 
iBerezinskv and DokuchaevI l|200l|) and the case (ii) we de- 
scribe below. In particular, in the case (ii) it must be taken 
into account a possibility of the SMS contraction prior to 
the formation of a central NS cluster. 

The NS dynamical friction time-scale for a mean NS 
inside the host SMS is decreasing during contraction of 
SMS as tdf oc _R^/2 ^^j^g contrary, the SMS evolution 
time-scale is increasing with SMS contraction as isMS oc 
R~^. Thus, the subsystem of NSs begins evolving faster 
than the host SMS after reaching the stage when t^t — 
tsMS- At this stage all NSs sink deep to a central part 
of SMS and form there a self-gravitating cluster. From 
relation tdf{R) — isMs(^) we find the corresponding SMS 
radius at the moment of NS cluster formation inside it: 

RsMS-{^ ^\ T ~4.6x IQi^Mg^/'cm, 2 

where Afg = Msms/(10'*Mo). Respectively a mean gas 
density 

3 MsMS „n inl7,,r2/5 -3 /o\ 

'T'SMS ^ 2~"p3 ~2.9xlO"M8' cm (3) 

a mean column density 

^SMS - nsMsmpRsMS - 2.2 x lO^Mg'^' g cm'^, (4) 

a virial velocity usms = (GMsMs/2i?sMs)^^^ ^ 
0.04Mg^^c, a SMS evolution time t^us — isMs(-RsMs) — 
7.3 X lO^Mg^^ yr and a surface black-body temperature 
TsMS — 3x 10''Mg^^° K. Note that the only free parameter 
in these relations is Msms- A self-gravitating NS cluster 



3 

with a radius i?NS Rsms is formed when a total mass 
of NSs inside this radius becomes of order of an ambi- 
ent gas mass, Mns — (47r/3)i?NS/°c, where pc is a central 
SMS density. From this relation the initial radius of NS 

1 /3 — 1 /3 

cluster is i?NS — /ns ^'^ Rsms and velocity dispersion 

j.l/3,1/3 , 

«NS = /ns ^SMS, respectively. 

4. Neutron star cluster evolution 

The formed NS cluster evolves much more faster than a 
host SMS and collapses finally into the MBH. A suitable 
approximation for the dynamical evolution of NS clus- 
ter in our case is the hornologou s "evaporation" model 
(e. g. ISnitzeil . Il987t ISaslawi Il987l) . This model describes 
the two-body interactions of stars in the cluster with an 
assumption that a fast star escapes (evaporates) from the 
cluster with a zero total energy. Accordingly, a total virial 
energy of the cluster remains constant E = —Nmv'^/2 — 
const during evaporation of fast stars. The velocity dis- 
persion in the evolving cluster is growing as w oc N~^^^ 
and cluster radius diminishes as i? oc N'^ with a dimin- 
ishing of number TV of remaining stars in the cluster. 
After reaching v ~ 0.3c, which is onset of the NS cluster 
global dynamical instabili ty, the remaining NS cluster col- 
lapses to the MBH (e . g. 'Zcl'd ovich and Podureti 119651: 
lOuinlan and ShaDirol[l987) ). The rate of NS evaporation 
from the cluster is 

iVev ^ aNt-\ (5) 

where a relaxation time t^ is given by and constant 
a ~ 10^'^ — 10^^ according to the numerical calculations 
in the Fok kcr-Plank approximation ('Spit zeilll987t ISaslawL 
Il987t ICohn. 1980). Integration of © together with the 
relation E = const results in the nondissipative evolution 
due to evaporation of fast stars: 

Nit) = A^NS (l - ^) ^ ' (6) 

where the cluster evolution time iov = Oovir, numerical 
constant Oov = (2/7)a^^ and t^ is the relaxation time at 
the moment of NS cluster formation t = 0. In numerical 
estimations we put = 10^02 with a2 ~ 1- 

The random collisions of NSs and their evaporation 
from the cluster occur at the same time. A rate of NS 
collisions in the cluster wit h the gravitational ra d iation 
losses taken into acco unt l|Quinlan and Shapirol Il987l: 
iDokuchaev et al.Lll998(l is 

iV..p^36V2(^)"''l-^, (7) 

where rg = 2GmNs/c^ is a gravitational radius of NS. 
According to © the collision rate {Tj) is gradually growing 
due to the evaporative dynamical evolution of cluster as 
A'cap v^^^'^ . The dissipative process of collisions is negli- 
gible for the dynamical evolution of a cluster in compari- 
son with the nondissipative evaporation until iVcap < iVev 
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This inequality is equivalent to the condition v < Ucap, 
where 



V7V3 flc 



(8) 



At u > t;cap the NS collisions becomes a major evolution 
process in cluster. The contraction of cluster is accelerated 
and it finally collapses to MBH. 

5. Fireballs and rarefied cavity inside the SMS 

We assume that each NS collision is accompanied by the 
generation of relativistic fireball with a total energy Eq = 
10^^i?52 erg, where £^52 ~ 1. The physics of fireballs is 
extensively elaborated in the recent years for the gamma- 
ray burst (GRB) models (for review see e. g. IPiranL I2OO0I 
and references therein). 

Although in our case at least the first fireball is born 
in SMS with large density nsMS given by Eq. the 
dynamics of the fireballs remains the same as in the stan- 
dard theory. Consider the first fireball assuming it to be 
spherically symmetric, which in fact is not a necessary 
assumption. 

The tremendous energy release Eq ^ 10^^ erg occurs 
within the small volume with radius Ri ~ 10^ cm and 
the radiation pressure pi ~ Eq/R^ by many orders of 
magnitude exceeds the external pressure p ~ nsMskT. 
The pressure in the fireball, which consists of relativistic 
7e"'"e~ gas with small admixture of the baryons, accel- 
erates fireball to large Lorentz-factor Tf. The relativistic 
shock, with Lorentz factor Tgh — f propagates ahead 
of the fireball pushing the SMS gas. The shock becomes 
non-relativistic at the Sedov length 



Is = (3So/4™sMsmffc')i/3 = 1.7 X 10^ 



(9) 



This value gives an order of magnitude estimate for cavity 
radius (see below). 

The second fireball propagates in the rarefied cavity 
produced by the first fireball, and it overtakes the first 
one at about Sedov length, when the first shock is non- 
relativistic. The shocks from the first and successive fire- 
balls merge into one, hence the following scenario is plau- 
sible. 

The numerous fireballs after successive NS collisions 
in a dense cluster provide the central energy source in 
our model. The power of this source is not limited by the 
Eddington luminosity in contrast to the accretion sources. 
The radius of NS cluster is very small in comparison with 
the host SMS radius -Rsms- The shocks produced by nu- 
merous fireballs dig out a rarefied cavity with a small 
radius in comparison with i?sMS- Consider the merging 
shock in more details. 

The early nonstationary stage of shock expansion can 
be described by the self-similar spherical shock solution 
for a central source with energy vary ing in time as £^ = 
At'' with A = const an d k = const ('Ostr iker fc McKeeL 
Il995l iDokuchaevL 120021) . The particular case of = 



corresponds to the Sedov i nstant shock solution (see e. g. 
iLandau and Lifshit j . llQS'^ . Meanwhile the nonstationary 
shock from multiple fireballs corresponds to a permanent 
energy injection into the shock, or injection shock with 
a central source of constant luminosity, k — 1, L — A, 
E = Lt. The radius of a self-similar expanding shock grows 
with time as 

i^ = i?(^) = /3(^^y^^(2+'=)/^ (10) 

where numerical constant /? — f3{'j, k) depends on the 
gas adiabat i c ind ex 7 and fc, e. g. /3(4/3,l) — 0.793 
llDokuchaevll2002|) . The velocity of shock expansion is 



dR 



{2 + k)R 
5t 



(11) 



The maximum radius of the expanding strong shock Rgh 
is obtained from the equality u{Rsh) = Cs.c by using H10|l 
and ((TT|l : 



Rs 



2+k 



5 Ci 



k\ ' 



A 



l/(3-fe) 



(12) 



Respectively the maximum time of a strong shock expan- 
sion is tsh = [(2 + fc)/5](i?sh/cs,c)- The expansion law for a 
fc = 1 injection shock corresponds to constant energy flux 
(or constant source luminosity) carried by the swept-out 
gas. This constant luminosity shock solution is applicable 
only to the early nonstationary stage of cavity formation 
at i < tsh- 

At the late stationary stage at t > tsh the boundary 
of the cavity is supported in a dynamic equilibrium by 
the relativistic wind from successive fireballs. The radius 
of a stationary cavity i?cav is determined from the energy 
fiux balance on its boundary. The central source power or 
luminosity is L — NcapEo, where Neap is from ((T)). After 
formation of cavity the energy generated in the central 
source is carried out by a relativistic wind. Just outside 
i?cav the energy fiux is carried by a hydrodynamic flow 
L = 4ttRI^^pv{w + w^/2) with a gas velocity v ~ Cs,c, 
where Cs,c and pc are respectively a sound velocity and 
gas density in the central part of SMS. The gas produces 
some work under expansion, so the energy flux L contains 
an enthalpy density w — e + p/p = €3/(7—1), where 
e — Cg 7(7(7 — 1)], p, p and 7 are respectively a gas inter- 
nal energy density, pressure, density and adiabatic index. 
From the energy fiux balance relation we determine a sta- 
tionary cavity radius for 7 — 4/3 



Rc 



N,: 



1/2 



lAnpcC^ 



(13) 



For the stationary cavity existence the time interval be- 
tween successive fireballs tc„^ 



i-cap — A^cap must bc less than 
that for cavity shrinking (or spreading) <cav — ^cav/<-s,c- 



This requirement, t 



cap 



with a help of and ifT^ 



can be written as condition v > Wcav, where 
z;cav 7.3 x 10"^ E~J/'''m!'^'''c, 



(14) 
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with numerical values Vcav — 4.3wns- The corresponding 
minimal radius of a stationary cavity i?min = -Rcav(wcav) 
is 



En 



1/3 



2.1 X lO^^E^s^^Afg cm. (15) 



It follows from itT^ that i?min is independent of /ns- At 
the moment of a stationary cavity formation inside SMS, 
V = Vca,v, we have respectively a radius of the NS cluster 
i?(i>eav) =^ 7.5 X 10^^f]_y^^El^/^''^Ml'/^^^ cm, a number 
of NSs in the cluster at the moment of cavity formation 
iV(t;cav) ^ 3.8 X lOV^y^^sa^'^'^^s''^'''^ a NS cluster 

, , . , . , / \ n o j-41/27 i:n49/135 , .rl73/225 

evolution time tcv(Wcav) — 0.3a2/_2 ^52 -'^'^s 

and central source luminosity inside a SMS i(wcav) — 

8.6 X 10^^i?52'^Afg^^ erg/s. After the formation of a sta- 
tionary cavity inside the SMS we have a hierarchy of radii, 
R{vcav) < Rmin ^ ^SMSi which justifics thc using of cen- 
tral values for a gas density pc and sound velocity Vs,c- 

The central source power and cavity radius are gradu- 
ally growing due to the evolutionary grows of NS collision 
rate. The rate of NS collisions and the power of the cen- 
tral source reach the maximum value at w ~ Wcap from 
||SJ|. The corresponding maximum central source power 

imax = A'"cap(Wcap)£'o is 



L„,ax ^ 2.8 X lO^^a-'/Vr'^ ^^52Mg-'/' erg s-\ (16) 

By substituting this luminosity in p3l) we find that the 
maximal cavity radius i?max — ^cav(i'cap) is 

D o o inl2 -9/4 7;il/2 , ,^-17/10 -n 

^max =i 3.8 X 10 flj /_2 -£^52 ^8 ^m. (17) 

Respectively at ?; = TJcap the NS cluster radius i?(fcap) — 
3.6 X 10^^02^^^/^2'^Mg^^ cm, number of NS in the cluster 
7V(wcap) ^ 2.7xlQ'^ay^fJ^Ml'^ and NS cluster evolution 



time tov(wcap) 



Zl.iaf"'fTMl"' days. 



Collisions of NSs during a cluster lifetime tcv at 



supply the total energy £^ 



7.6 X Kf'^al'^ fJ^EoMl'^ erg into the cavity. This en- 
ergy is far less than a SMS binding energy ii^sMS — 
GM|Ms/2i?sMS - 2.9 X IQ^'^mI'^ erg, and so the shocks 
after numerous fireballs do not influence the SMS internal 
structure. 



6. Production and detection of high-energy 
neutrinos 

We shall study here the acceleration of protons by rela- 
tivistic fireballs produced by NS collisions and the high- 
energy neutrino radiation produced in p7 and pp colli- 
sions. We assume a standard GRB fireball which prop- 
agates at the baryon-dominated stage with the Lorentz- 
factor r = 10^r2, r2 ^ 1 and carries a baryonic mass 
Mq = EqIV^ ~ 5.6 X lO-5£;52r2"^M0. 

The acceleration occurs due to internal shocks inside 
the fireball, iWaxman. (,1995^1 , and due to external shocks. 



IVietril l)l995|) . We shall consider first the internal shock ac- 
celeration and high-en ergy neutrino production following 
IWaxmanI l200lh . This scenario we consider as the 

basic one for a prediction of neutrino flux in our model. 

Th e standard explanation of the observed GRBs. lPiranI 
l)20fl(]h . is given under an assumption that fireball con- 



sists of the mildly relativistic sub-shells. The inner sub- 
shells are moving faster than the outer ones, and their 
collisions produce mildly relativistic shocks. The protons 
are accelerated by these shocks with the standard spec- 
trum cx i?^^ and wit h a maxi r num a cceleration energy 
i^max - 1 X 10^1 eV, IWaxmanI l|l995l) . The specific fea- 
ture of our case is a relatively short time of fireball ex- 
istence in the cavity, t ^ i?cav/c. This time, however, is 
long enough for acceleration of protons up to the maxi- 
mum energy E'max ^ 10^^ eV. In the rest-frame K! the 
acceleration time is t'„ ^ r r / c. where is the Larmor ra- 
dius. Following to lPiranI l)200nl) . an equipartition magnetic 
field B' in the relativistic wind can be estimated as 



n'2 



2Ci 



(18) 



where unprimed symbols correspond to quantities in the 
laboratory frame, io'^ is photon energy density, ~ 1 is 
an equipartition parameter and L is a fireball luminosity. 
Using L-y ~ 10^^ erg/s and r ~ i?cav we obtain the acceler- 
ation time as ia ~ 50 s for particles with E'max ~ 10^^ eV. 
With an equipartition magnetic field (|18|1 the maximum 
acceleration e nergy is E ~_eB' i?cav ~ 10^^ eV in agree- 



ment with the IWaxmaiil ( model 

High-energy neutr inos are produced in P7 collisi ons in- 
side the fireball by the lWaxman and BahcalJ (|1997^ mech- 
anism. Due to threshold of pion production in p^ collisions 
the neutrino flux has a low-energy suppression. The low- 
energy part of neutrino spectrum is produced in pp colli- 
sions of accelerated protons with a gas inside the cavity 
and in the SMS interiors. Reflecting from the fireball, the 
cavity gas penetrates into the fireball at the Larmor radius 
distance d ~ r^, for a particle with energy E' ~ Fmp. In 
this region accelerated protons interact with an accumu- 
lated gas, producing pions and neutrinos. However, below 
we shall consider more efficient mechanisms for production 
of low-energy neutrinos. 

We assume, as it is usually used in the GRB cal- 
culations, the equipartition of total energy distribution 
in the form of UHE protons W^, neutrinos Wy, GRB 
photons VFgrb and the energy transferred to the shock 
in SMS. Then for neutrino luminosity during the max- 
imum activity of the considered central source we have 
Ly ^ (l/4)L,„ax with Lmax giveu by (fT?)|l . Though neu- 
trino spectrum is presented by two components, for sim- 
plicity of further estimates we assume the universal spec- 
trum oc 1/ E"^ . Normalizing this spectrum to ij^, we obtain 
for the neutrino generation rate 



1 



0. = 7 



4 £2 \n{Ey^^,^/E 



(19) 
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Taking into account the neutrino oscillation, one obtains 
^'p + j^fi neutrino flux at a distance r from the source 
J^..+..(i^) = (l/3)0.(i?)/W. 

We shall discuss here only the detection of TeV neutri- 
nos from the considered hidden source by a future 1 km^ 
ice or underwater detector. The fluxes of UHE neutrinos 
up 10^^ eV are most probably also detectable, but they 
need the knowledge of the details of neutrino detection. 
This task is beyond the scope of this paper. 

Crossing the Earth, neutrinos with energies E > 0.1 — 
1 TeV produce deep underground a n equilibrium flux of 
HE muons, which can be calculated ijBerezinskvL Il99(i|) as 

1 



F^{> E) 



(20) 



where ctq ^ 1 
section, Na = 



X 10^34 ^j^. 
6 X 10^3 



^ is the normalization cross- 



is the Avogadro number, 6^ 



4 X 10~^ cm^/g is the rate of muon energy losses, ( = 
In i^max/'E'inin, and Yii.(E) is the integral niuon moment 
of ly/ .N interaction ((Berezinskv et all llflflOt iBerezinskvl 
Il99n(l . The most effective energy of muon detection is 
Efj, > 0.1 — 1 TeV, because these muons cross the large 
part of 1 km^ detector, and their directions can be reli- 
ably measured. The number of muons N{v^) — Ff^Stcv 
with > \ TeV detected during the burst duration 
^ov = 31.3 days by the underwater detector with an ef- 
fective area S at distance r from the source is given by 



N{v^) ~ 56 



in 



S 



2.8 X 1049 erg s- 
The detected muons with E,, 



(21) 



1 km" Mpc, 

j^>l TeV show the direction 
to the source with an accuracy 2° — 3°. 

The other acceleration mechanism due to the external 
ultra- relativistic shocks, iViet ri ( 1995.) , provides mostly pp 
neutrino flux with energies below 1 TeV. Consider the 
shock which is produced when a fireball hits the inner 
surface of SMS. The shock propagates in ultra-relativistic 
regime at distance of order of the Sedov length, before it 
overtakes the slow main shock and merges with it. The 
gas density in this region rises from that of the cavity 
Tie ^ 4x 10^"'^ cm~^ to maximum density rimax — 54.2nsMSj 
with nsMS — 3 X 10^'' cm~'^ given by The high- 
energy proton initiates a nuclear cascade due to pion pro- 
duction. The energy dissipates from the cascade through 
cascade-pion decays. The decays of neutral pions initi- 
ate the electro-magnetic cascade and the released energy 
contributes to the shock and heating of SMS. The col- 
umn density of SMS, Xqms - 2.3 x lO^Mg^^ g cm"^, 
is large enough to absorb all produced particles except 
the secondary neutrinos. Neutrinos are produced in de- 
cays of the charged pions. The critical Lorcntz factor of 
the decaying charged pion is determined by condition ^ 
[aT^tqUTT^c)"^ ~ 5 X 10^, where cTttA' — 25 mb is ttN cross- 
section and — 2.6 x 10~* s is the charged pion lifetime, 
and for the mean gas density we used n ~ 1 x lO^'* cm^"^, 
larger than nsMS and smaller than Jimax- As a result we 
obtain for the maximum energy of the secondary neutrinos 
i^.,max ~ 1 TeV. 



Acceleration by ultra-relativistic shock at the inner 
boundary of SMS proceeds in a specific way because of 
the large density of matter there and large magnetic field. 
Consider one of the fireballs hitting the inner boundary 
of SMS. Propagating through the dense gas the shock, 
running ahead of the fireball, accelerates the protons. A 
proton from the SSM interior (upstream region), which 
makes its first downstream - upstream - downstream (d- 
u-d) cycle, i.e. crossing the shock and scattering upstream 
back to downstream, increases its energy by a factor 
Fgj^. In the co nsequent d-u-d cycles the ene rgy gain is 
only about 2, 'Gallant and Achterberd l)l999f) . (see also 
iLemoine and Revenu . (,2005j) for a recent work) . The max- 
imum acceleration energy is determined by the number of 
completed d-u-d cycles and is given by ta ^ iint condi- 
tion, where ta and imt are acceleration and interaction 
time, respectively. They are g iven by ta ^ E/ecBTsh, 
lOallant and Achterberd l)l999(l , and iint ^ 1/ctTppnsMS- 
Acceleration time is short because of very large magnetic 
field accumulated at the inner boundary of the shock. 
The magnetic field in a fireball, Boq, given by Eq. H18|l 
reaches 4 x 10^ G, when it hits the SMS. The each succes- 
sive fireball adds the magnetic flux, and though not be- 
ing summed up coherently, magnetic field increases. For 
_B ^ 4 X 10^ G the maximum energy i?max ~ 1 x 10^'' eV. 
Fortunately, this energy, being unreliably estimated, is not 
important for our calculations of neutrino flux, because all 
high-energy protons loose energy by producing nuclear- 
electromagnetic cascade. 

The proton with an energy Ep transfers fraction of 
energy f{Ep) to neutrinos through decays of cascading 
charged pions, T,Ey = f{Ep)Ep. The fraction f=l/2 at 
Ep < 10^^ eV and correspondingly / ~ 0.05 at Ep ~ 
10^^ eV, because in the cascading process all neutral pi- 
ons decay. For the total energies of protons and neutrinos 
Wi, — fWp, where / is the value averaged over spectrum 
of accelerated protons. We will assume / ~ 0.1. It is easy 
to see that neutrino signal in the detector is determined by 
W^. It occurs because the cross-section at the energies 
of interest is a{E) = ctq^^, with ctq = 0.8x10"^* cm^/GeV. 
The total number of neutrino events in 1 km'^ detector is 
given by 



N^{E)EdE 



(22) 



where the number of target nucleons is Nn — 6 x 10^^ 
and ~ 0.1 Wp. For r = 1 Gpc and Wp = (l/4)iytot = 
1.9 X 10^^ erg, it results in Atot = 47 of low-energy {E^, < 
1 TeV) neutrino events in 1 km'^ detector. The neutrinos 
produced in the SMS by high-energy protons, accelerated 
in the inner shocks, double this number of neutrino events. 

On may wonder about the influence of adiabatic en- 
ergy losses of accelerated protons on neutrino fluxes, esti- 
mated above. The role of adiabatic energy losses has been 
critically dis cussed by iRachen and MeszarosI l)l998l) (see 
also reply bv lWaxmanI l|2004l) . For the calculations above, 
the adiabatic energy losses are irrelevant, because in all 
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cases neutrinos are produced in situ with the acceleration 
of particles. 

How many these hidden sources are expected? A fur- 
ther (after neutrino burst) evolution of SMS with a mass 
M ^ IO^Mq inevitably results in formation of a mas- 
sive black hole and the AGN activity. Considering the 
described hidden source as a "pre-AGN" stage one can 
estimate an expected number of hidden sources inside the 
cosmological horizon. Since a hidden source and AGN are 
considered as the different stages of the same evolving 
galactic nucleus, the probability to have a hidden source 
in the set of AGN is given by the ratio of durations of the 
corresponding evolution phases, tev/^AGN, where tAGN is 
a duration of AGN phase. Then the number of hidden 
sources within the horizon is given by 

A'^HS ^ g7r(3rfo) "-AGN U/tAGN , (23) 

where (47r/3)(3cto)^ is the cosmological volume inside the 
horizon ctg and jt-agn is the number density of AGN. The 
estimates for wagn and ^agn taken for different popula- 
tions of AGN results in iVns ~ 1 - 10. 

Can the hidden sources be observed in electromag- 
netic radiations? X-ray and gamma-ray radiations are 
fully absorbed due to tremendous thickness of SMS, 
-'^SMS '^-^ 2 X 10^ gcni^^. However the thermaized radi- 
ation can propagate diffusively and reach the the outer 
surface of SMS. This process h as been considered in 
iBerezinskv and DokuchaevI l)200lj) for the case of a mas- 
sive shell. In case of SMS, the diffusion is much slower and 
the thermalized radiation reach the outer surface much 
later, and thus will be emitted with a lower luminos- 
ity. The diffusion coefficient D cl^, where a diffusion 
length is Id = 1/((Tt'^sms)i f^nd ctt is the Thompson cross- 
section. A mean time of the radiation diffusion through 
SMS is very long, ta RIms/^ ~ 4 x lO^Mg^^ yrs, in 
comparison with a lifetime of the NS cluster inside the 
SMS. Since the dispersion of the diffusion time distribu- 
tion a ^ td, the surface luminosity of the central source 
is Lbb ~ £^tot/o- ~ L^B^^U^r/td 5 X 10^^ erg/s and hence 
is negligible in comparison with the Eddington luminosity 
of the host SMS. 

7. Conclusions 

In this paper we studied the natural scenario for massive 
black hole formation in the process of dynamical evolu- 
tion of a compact stellar cluster. Our contribution to this 
scenario is based on the observation that an intermediate 
SMS can be formed with a subsystem of neutron stars sub- 
merged in it. The evolution of neutron star cluster results 
in its contraction and merging of neutron stars at their 
collisions. Both these processes have been studied before 
and rather well known. The merging of neutron stars can 
produce quasi-spherical fireballs, where protons are accel- 
erated to high energies, and neutrinos are produced. 

Described in a more detailed way our scenario looks 
as follows. The dynamical evolution of a central stellar 



cluster in the galactic nucleus results in the collision de- 
struction of normal stars and formation o f massive gas 
envelope l|Berezinskv and nokucha,ev[l2nnih or SMS (this 
paper). The neutron stars and stellar mass black holes 
survive through this stage and form self-gravitating sub- 
system deep inside the SMS. The lifetime of this compact 
neutron star cluster is very short. The multiple fireballs 
after neutron star collisions in this cluster dig out a rar- 
efied cavity inside the SMS. The protons are accelerated 
in the fireballs and in the cavity and produce neutrinos 
in collisions with an ambient photons and gas. The ac- 
tive phase of high luminosity of the central source has a 
duration t ~ 0.1 — 1 yr. During the neutrino burst the 
energy of the same order as a total neutrino energy 
is released also in the form of high-energy electrons and 
photons. All high-energy particles, except neutrinos, are 
completely absorbed in the SMS, which column density 
can reach X ^ 10^ g/cm^. Thus the considered object 
is a very powerful hidden neutrino source. It is demon- 
strated that that shocks after numerous fireballs are ab- 
sorbed inside the SMS and cannot destroy it. The system 
remains gravitationally bound and gravitational collapse 
of the compact neutron star cluster produces a seed mas- 
sive black hole inside a host SMS. Long after the neutrino 
burst, the described hidden neutrino source is to be seen as 
a bright AGN. As precursors of most powerful AGN, these 
hidden sources are expected to be at the same redshifts 
as AGNs. The number of the hidden neutrino sources can 
be estimated as ^ 1 — 10 in the visible universe. 

The hidden neutrino source is detectable by future 
1 km'^ detector with a number of TeV muons up to 
50 — 100 per source. The predicted number of low-energy 
{El, < 1 TeV) neutrino interactions is ~ 100. The high- 
est energy neutrinos with E > 10^^ eV are probably also 
detectable. 
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